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Abstract 

We discuss a possible strategy for studies of a particular scenario of New 
Physics (NP) at LHC. The NP is taken to be a). C/-spin symmetric, i.e. it does 
not distinguish d and s quarks; b). it makes no contribution to the tree processes, 
but contributes differently to penguin and box diagrams and c) . it does not spoil 
the unitarity of CKM matrix. Our analysis is based on comparison of particular 
CKM matrix elements, which can be obtained from the processes dominated by 
diagrams of different topology. We argue that the standard formalism of the 
overall unitarity triangle fit is not suitable for studies of NP of this kind. We 
also stress the utmost importance of lattice computations of some particular set 
of hadronic inputs relevant for NP searches. 



1 Introduction 



The main interest of the present day flavor physics is focused towards search- 
ing for possible signals of New Physics (NP) - the effects which are not taken 
into account by the Standard Model (SM). These still hypothetical effects can 
be roughly divided into two groups. The first, "quantitative" one consists of 
effects which are present in the SM but whose concrete SM prediction deviates 
from actual experimental results. A well known example is given by rare de- 
cays strongly suppressed in the SM but expected to be enhanced in some NP 
scenarios. Another, "qualitative" group, is formed by the effects which are not 
present in the SM at all, like possible observation of nonconservation of any 
charge (baryon, electric etc) strictly conserved in the SM. At the moment there 
are no clear indications on possible NP effects of either kind. The strong hope 
however is that the situation will change in the nearest future with the run of 
LHC. Of prime importance in flavor physics is an analysis of the CKM mixing 
matrix. The commonly accepted parametrization-independent language used to 
discuss the rich physics encoded in CKM matrix is formalism of the unitarity 
triangle (UT). For introduction into the subject and all details the reader is 
refered to the materials presented on [U [2l [3l |4] and to excellent recent reviews 

mm- 

The issue of NP search in flavor physics context is certainly much broader 
than the mere check of CKM matrix unitarity, however precise it can be. Of 
course, any inconsistencies in the UT construction will undoubtedly indicate the 
presence of physics beyond SM. The opposite is far from being true - there are 
many reasonable NP scenarios which are well compatible with perfect unitarity 
of CKM matrix. 

There are different possible strategies to study the CKM matrix. The most 
popular one, adopted in particular by UTfit and CKMfitter groups [3 12 is to 
use all available experimental data to overconstrain the triangle. Besides general 
importance of this activity the hope is that the procedure will exhibit some 
inconsistencies signaling NP effects. Up to now there is an overall agreement of 
all constraints (see recent talks [Illl]). 

However, this approach also has some disadvantages. In our view, the most 
important one is the fact that the set of constraints in use is not fitted to this 
or that particular NP scenario. On the other hand, the relevance of this or that 
observable from the point of view of its possible NP content strongly depends 
on what kind of NP we discuss. Let us explain this point taking as a typical ex- 
ample AMg/ AMij ratio. For all scenarios where NP couples identically to s and 
d quarks ([/-spin symmetric NP) this ratio is not sensitive to NP contributions, 
since in this case short-distance functions, even if modified with respect to the 
SM predictions for each AM^, AMg exactly cancel in the ratio. This pattern is 
typical for, e.g. constrained minimal flavor violation (CMFV) NP models (see 
review of MFV models in [H]). As a result this quantity informs us about ratio 
of couplings of t-quark to d and s-quarks and also about long-distance SU (3) 
breaking effects in QCD (see expression pUj) below), but brings no information 
about correctness of the short-distance SM calculation of AMd or AAIg sepa- 
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rately. And it is precisely the latter short-distance piece we are interested in 
most of all if we are looking for deviations from the SM at small distances. On 
the other hand, there are NP scenarios such as, e.g. MSSM at large tan/3 (see 
[S] and references therein) and next-to-minimal flavor violation [10] where this 
is not the case and the ratio under study is sensitive to NP. Moreover, it is very 
natural to expect (and this is our general attitude in the present paper) that 
NP contributes differently to the processes of different topology (i.e. tree and 
penguin, penguin and box etc). Obviously, this effect can be lost in compari- 
son of observables of the same topological type. In view of that an alternative 
way has been proposed (see, e.g. [121 US] and also [TJ [51 [SI [TH [TS] and refer- 
ences therein). Generally speaking, it corresponds to construction of a few a 
priori not coinciding unitarity triangles, each extracted from branching ratios 
and asymmetries for processes of some particular kind. In this case any mis- 
match between these UT's, e.g the so called "reference UT" [13] and "universal 
UT" (see recent discussion in [9]) would be a clear signal of NP, and, moreover, 
one could in principle identify the place (EW penguin sector is among the most 
promising ones) where it has come from. 

Adopting the basic idea of the latter strategy we address the following prob- 
lem. Let us assume the still hypothetical NP is, in the spirit of next-to- minimal 
flavor violation scenario: a), [/-spin symmetric; b). does not contribute to the 
tree processes and c). does not spoil the unitarity of CKM matrix (i.e. we work 
in the spirit of next-to-minimal flavor violation scenario). How can we see NP 
from global UT fits and what observables are the most sensitive to NP effects 
in this particular case? 

To answer this question, we analyze theoretical and experimental (having in 
mind mostly the LHCb experiment) perspectives for studies of some CKM ma- 
trix parameters which can be extracted from the processes of different topology 
and can be sensitive to NP of the discussed type. It is worth noticing that the 
mismatch between sin 2/3 values from B — > J/ipKs and from B — > (f>Ks modes 
widely discussed in recent literature (see, e.g. [5] [TB] and references therein) 
represents exactly a kind of effects we are interested in. We will also stress the 
urgent need for new refined lattice data on hadronic input parameters in order 
to determine the product |VtsV(J|. 

The paper is organized as follows. The section 2 is devoted to brief overview 
of the existing strategies for CKM matrix analysis, while our procedure and 
results are presented in the section 3 and conclusion in the section 4. 

2 Overview of the standard strategy 

In general one can choose different sets of independent parameters which enter 
the basic unitarity relatior0 

Vudv:, + VM + Vuv;, = (1) 

^As well as five other unitarity triangles 
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It is worth noticing that the term "independent" is usually used in the litera- 
ture in mere algebraic sense, i.e. one assumes no relations between CKM matrix 
elements other than those following from the unitarity constraints. This assump- 
tion may be wrong if some more fundamental underlying structure behind CKM 
matrix does exist. A common choice for one of the parameters is 

_ r (0.2265 ±0.0020) P , 

si2 - A- - I (Q 2258±0.0014) [2] 

This quantity can be determined with very good accuracy from the decay mode 
K ttIv with the latter being dominated by tree level process. The main 
source of error here is the poor knowledge of the corresponding formfactor /+(0), 
namely, according to [H] (5|Ks|/_^(o) = ±0.0018, S\Vus\exp = ±0.0005. 
The interior angles of the triangle ^ are conventionally labeled as 

The Cabbibo-suppressed angle x important for Bg — Bg oscillations is defined 

by 

and is also of interest. Let us briefly remind the strategy for 7. The cleanest 
way to extract it is from the interference of the b cus and b — > cus transitions 
(the so called "triangle" approach). Practically, this corresponds to the study 
of B~ — > K^D^ and B~ K^D^ modes with the subsequent analysis of the 
common final states for D and D mesons decays. One considers CP-eigenstates 
as final states for D,D mesons decays (GLW approach [18]) or combines ob- 
servables from different modes {B K*D, B KD*, B KD, B K*D*) 
(ADS approach [19]) to overconstrain the systcmH Notice that the interfering 
diagrams are the tree onesH 

The combined results for 7 presented in [3] obtained by Dalitz plot analysis 
|21| are given by 

7 = 67° ± 28° ± 13° ± 11° [BaBar] ; 7 = 53° ± 18° ± 3° ± 9° [Belle] (5) 

where the errors are statistical, systematic and the error resulting from the 
choice of Z? - decay model. For the discussion of the situation with 7 determi- 
nation in LHCb the reader is referred to [22] . 

Using various methods the overall uncertainty in 7 at LHCb is expected to 
be as small as 5° in 2fb^^ of running and will eventually reach the level of 1° 
with increase of statistics. 



^The same strategy can be applied to the case of Be mesons, where it has some theoretical 
advantages I20| : however the experimental statistics becomes the main obstacle there. 

•^This mode is not completely NP-safe since in principle the latter can enter through D^—D^ 
mixing. 
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With the standard assignment for the elements of CKM matrix (see, e.g. 



E) 

si2 = A ; S23 AA^ ; S13 exp(-i(5i3) = AX^{p - ir]) 
to define the apex of the unitarity triangle 



1-U 



l-^A^ 



(6) 
(7) 



one needs to know at least two independent quantities out of two sides 



Rb 



\Vudv:,\ 



VcdV^ 



Rt 



\VtdV,l\ 



cfcl 



'cb\ 



and three angles a,/3,7 where the latter are defined by ([3]). In particular, 
the authors of [15] analyzed all ten possible strategies, distinguished by the 
mentioned choice of two independent parameters out of five from the point of 
view of their efficiency in the determination of UT. For example, our geometrical 
intuition tells us that it is easier to construct general non-squashed triangle 
taking as inputs one of its angles and adjacent side (because the variations in 
these parameters are approximately orthogonal) than taking the same angle and 
the opposite side (because the variation in these parameters are approximately 
parallel) . Numerical simulation done in }15j fully supports this intuition, giving 
the highest priori tj0 to the strategies based on combined use of either (7,/?) or 

This result is particularly encouraging because the quantities Rb and 7 define 
the so called reference UT [121 US] • The latter is built from the observables that 
are expected to be unaffected by NP, since their dominant contributions come 
from tree level processes. Then assuming unitarity of CKM matrix one can 
compute from ([1]) reference values for 



Rt 



1 + i?? - 2i?f,cos7 ; cot/3 = 



1 — Rb cos 7 
Rb sin 7 



(8) 



and compare them with the ones obtained by direct measurements in the pro- 
cesses involving loop graphs. Any difference could be a hint for a NP signal (see 
recent quantitative discussion of this issue in [9]). 

The elements of CKM matrix which enter the definition of Rb (up to terms 
0{X*) 



Rb 



1 Kb 



A \Vcb\ 



are known from semileptonic B-decays. The recent inclusive update is given by 
[1] \Vcb\ = (41.79 ±0.63) ■ 10-3. 

Experimental determination of |Kib|mc! suffers from uncertainties, intro- 
duced by specific cuts one has to apply in order to get rid of & — > c background. 



*In other words, to compute (p, rj) with a given precision pair (7, /3) may be known with 
lower accuracy than, e.g. pair (ij^,/?). 
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As for |V^f,|exc! the main source of error is lattice uncertainty in calculations of 
B ^ TT, p form- factors. Up to date results are given by [3] as 

|K6|»d = (4.4±0.3)-10-=^ 

and 

|K6|e.c/ = (3.8 ±0.6) -10-3 

At the moment the perspectives to increase the accuracy in experimental 
determination of Rb up to a few percent level are unclear. As can be seen from 
Fig.l and Fig. 2, the errors in 7 and Rb play a very different role in fixing the 
angle /3 with some given precision, which is a simple consequence of the fact 
that the angle a is close to 90° and the triangle is almost rectangular. The 
present accuracy in /3 extracted from the " golden mode" B ^ J/ ijjKs is better 
than ±2°, the current world average for sin 2/3 from tree level decays provided 
by [3] is (see recent talk [53] and references therein): 

sin 2/3 = (0.674 ±0.026) (9) 

The corresponding penguin contribution to f3 is Cabibbo-suppressed (see, e.g. 
[TB]V As shown in Fig.l an uncertainty window of ^ 3° for P corresponds to 
an uncertainty window of ^ (24 ± 5)° for 7 and therefore the precise data (O 
does not constrain 7 via ([5]) strongly enough to make the comparison discussed 
above meaningful. On the other hand, since both (3 from ^ and 7 from ([5]) 
are determined from the processes dominated by tree level decays, we do not 
expect to see violation of the second expression from ([5]) with these values of the 
angles. Anyway, the experimental uncertainty in 7 and hadronic uncertainties 
in Rb make ([5]) not valuable. 

Let us briefly discuss the side Rt- The are two ways of extracting Rt by 
means of relations not affected by NP contributions in some scenarios, notably 
CMFV. These are the computation of Rt from the first expression in ([8]) and 
the computation from the ratio AMd/ AMg where, again, short distance contri- 
butions to the box diagrams are canceled^ Concerning the former algorithm, 
because of the same geometrical reasons (angle a close to 90°) Rt is sensitive 
to the uncertainty in the angle 7 only (see Fig. 3). Thus, precise knowledge of 7 
will constrain Rt effectively. In the latter approach one obtains the ratio 



\Vt,\ "V"iB. AM 
with the nonperturbative parametejfl 



uh^^ (10) 



e = ^ (11) 

iJB. 



^As has been already mentioned, in MSSM at largo tan /3 the quantity Rt is sensitive to 
the different Higgs couplings to d and s quarks. 
^5 = 1 in case of exact flavor 5(7(3). 
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The typical error of current lattice simulations of ^ is estimated as 6% (see 
and recent analysis in [25]). Since up to 0{\'^) 



then having at our disposal recent CDF results [Hj (see (HH)) we can straight- 
forwardly extract for the mean value Rt = 0.92 with the uncertainty dominated 

byC- 

Let us summarize this part. Suppose we would be able to measure Rb and 7 
with some very high precision. This defines the position of the UT apex which 
is universal as soon as NP does not contribute to tree processes Rh and 7 have 
been extracted from. Let us also assume that we get Rt from AM^/AMs and 
/3 from B J/ipKs, and these observables perfectly agree with Rb and 7 via 
([U (i.e. the UT apex defined from Rt and /3 coincides with the one founc(l| from 
Rb and 7). Does this fact mean dramatic shrinking of NP parameter space? 
Not at all: for NP scenarios with [/-spin invariance and without sizeable NP 
mixing effectii this coincidence is trivial and brings no any information about 
the parameter space. One can tell that UT is simply too rough tool to see NP 
of this kind. In other words, the precise knowledge of ^ is important in this case 
to calibrate the lattice, but not to find the NP. 



3 Direct comparison of CKM matrix elements 
from different processes 

In what follows we are going to explore a complementary strategy whose essence 
is the comparison of values of CKM matrix elements, obtained from processes 
with dominant contributions coming from diagrams of essentially different topol- 
ogy. Again one can consider angles and sides in this respect. We are interested 
in observables, corresponding to the processes whose dominant contributions 
come from topologically different diagrams, namely: 

• radiative penguin in decay modes B —> K*"f, Bg — > 07 and B — > [p,uj)^, 
Bs K*j, for s and d quarks, respectively 

• oscillations of neutral B^ and Bg mesons with the dominant contribution 
given by box diagram, resulting in the mass shifts AMg, AMd 

• tree and strong penguin interference in B decays into 2-body final states 
made of light hadrons tt, K , p and mixing relevant for the angle a deter- 
mination 

^Of course, any other pair can actually be used, see discussion above. 

*NP physics contributions in the box diagrams could in principle affect both 7 and /3 via 
^0 ~ Dq and Bq — Bq mixings. 
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For the first and the second mode the object of our interest is the product 
|V^^Vts|. For the third mode we confine our attention to the angles a, (3 and %. 

The reference values for these quantities are defined from tree level processes, 
since we adopt the usual assumption that they are free from NP pollution. One 
can make use of the "tree level definition" for [Vj^l^sl ^ |V'j^Vfs|tree, which up 
to terms ©(A"*) reads 



\VtsVtl\tree ^ \Vcb\ 



2 



l-y(l-2i?6COS7) 



(12) 



We have already discussed the corresponding numerical values and their uncer- 
tainties. Plugging them in, we get 

|^t.^t*fc|tree = (41.3 ±0.8) -10-3 (13) 
As for the angle a, its reference value is given by 

atree = TT - - 7 (14) 

where the extraction of P and 7 from tree level processes is described above. 
As for the angle %, there are no experimental constraints on it at the moment. 
The SM prediction is |x| ~ 0.02 0.04. 

3.1 Analysis of 

Let us start with the analysis of iVj^Vtsj. Values of these elements of the CKM 
matrix must exactly coincide in the SM, regardless of the way they are extracted. 
On the other hand, lack of such coincidence will be a definite signal of NP, 
contributing differently to these different types of processes. Qualitatively, one 
can consider ratios of the following kind 

(1) ^ \VtlVtq\AM^ (2) ^ \Vtb^tq\B^V-, A3) ^ \Vtb^tq\tree ,^r. 

where q = d, s and V stands for K* , (f>, p, lo. Thus we have three ways to extract 
the product |V(JVtg| of CKM matrix elements: via expression ^TE\i from the 
process dominated by the radiative penguin diagram, via expression (jl9p from 
the process dominated by the box diagram, and via (|12p from the reference tree 
level processes. It is obvious that by construction one has 

C • C • C - 1 (16) 



C^ = ei = C^l = ^ (17) 



In the SM however much more restricted condition has to be fulfilled: 

It is convenient to present a set of three numbers {Cq^v , Cq% as a single 
point on the ternary coordinate system with log C^'y as an (algebraic) distance 
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from the i-th axis. Then the SM case corresponds to the only pomt on this 
diagram - its origin, while any deviation from it is a hint to NP. 

The analysis of ratios of the mass shifts AA/^/AMs and branchings Br{B — s- 
pj)/Br{B — > K*j) widely discussed in the recent literature [531 H?! US] deals 
in our language not directly with the quantities Cqv^ but with their ratios like 

Cs^K*/C(il- '^be important advantage of these ratios is the improved accuracy 
of their theoretical determination, especially from the point of view of hadronic 
uncertainties. However the price to pay is high - the short distance factors which 
could contain contributions of NP are canceled in these ratios. In logarithmic 
coordinates it corresponds to a parallel translation, which could miss a consid- 
erable piece of NP, which is clearly seen from the analysis of [25]. In short, 

= Cd]l = 1 imphes Cl^A:*/Cd!p = 1' but not vice versa. 

Generally speaking, it is meaningless to look for (short-distance) deviations 
from the SM predictions if one has no quantitative knowledge what the latter 
actually are. Therefore as soon as we are discussing absolute values of mass 
shifts, widths etc, these short distance parameters should be determined and 
not just canceled in the ratios. Corresponding loss in an accuracy for hadronic 
contributions is perhaps inevitable. Anyway we are stressing that one has to 
deal with this " less accurate" low energy hadronic inputs if one tries to capture 
the short-distance effects of NP. For example, it is meaningless, in our view, 
to consider soft quantities as free parameters to fit observable branching ratios. 
Any possible NP induced difference between, e.g. the SM prediction for Br{B 
V^) and actual experimental result would be just hidden inside such "extracted 
from experiment" ''(0)1, which is clearly unacceptable. Simply speaking, to 
discuss deviations from the SM prediction we have first to know the latter □ 

In principle, one can discuss five expressions of the kind (fT6|) . corresponding 
to the following choices for (g, F): [s,K*), {s,(j)), (d, w), {d,p), {d,K*). How- 
ever all these channels have universal short-distance structure, while the long- 
distance contributions are related to each other by SU{3) flavor arguments. The 
optimal strategy therefore seems to choose just one particular case, which we 
take to be (s, K*) in the rest of the paper. The results for the other ones could 
provide important cross-checks (like, e.g. |Vtd|/|Vts| ratio), but presumably no 
new information about a NP content of ^TE\i . 

We are using the standard SM expressions for the decay rate for B K*^ 
and the mass difference AM^. The former can be written as [28l l29l [30| : 

ns - if*7) ^ ^^3,^!"' (1 - rf\a,{pWpm'\VtsV:,\^ (18) 

In the above expression r = m^*/m^, to;, stands for the pole mass of 6-quark, 
aj{fj) = Cy°^ -I- A(-^'(/x) is an absolute value of the corresponding short-distance 
function including Wilson coefficient Cy"' , hard scattering contributions and 

^As is correctly pointed out in [2] "a model-independent UT analysis beyond the SM 
cannot be carried out without some a priori theoretical knowledge of the relevant hadronic 
parameters." 
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annihilation corrections. The detailed computation of this function at the next- 
to-leading order can be found in the cited papers. Notice that we omit terms 



of the order of nig /ml 



The factor |^^^*^ 



differs from the corresponding 
form-factor T^^^ (0) by 0{as) corrections; according to [ST numerically one 
(0)1 « 0.93- rf^^*(o). 
The expression for /S.Ms reads as follows: 



has 



G 



AM, = ^i^B[M'^Fu]mBABBjl^\VtsVt 



* |2 
tb\ 



(19) 



where rjB is calculable short-distance QCD factor, while the mt /Mw-dependent 
factor FttM'^r has come from calculation of the box diagram ([32], see also 

[MlEl])- 

According to our strategy we invert the expressions (fT5)) and (fTO|) to the 
following form: 



WtsV, 



and 



WtsV, 



47r^ 

|a7(M)l 



2V{B K*j) 



G\am 



1(1 



rf 



kr^(o)i 



rub ) 



6AMs 



tsVtblAM, 



\/VbFu 



fBA Be. 



(20) 



(21) 



The structure of the above expressions is clear. The first factors in the 
r.h.s. are the short-distance SM contribution, which have to be calculated an- 
alytically. These are just numbers of order 1 and it is assumed that we have 
reliable theoretical control of this part. The typical accuracy of these factors 
is better than 5%. The second factors (the square roots) are composed from 
experimentally measurable quantities. The error in these factors is dominantly 
experimental and is currently at the 5% level for ((20l) and 1-2 % level for ((2T|) . 
The third factors (in the square brackets) encode information about soft QCD 
contributions (and related problem of b quark pole mass rub) for which we have 
no systematic approach of studying. The main hope here is focussed on the 
lattice simulations F°l The uncertainty of currently available data can be conser- 
vatively estimated as 10-20 %. The use of ((20)) . (|21l) as probes for NP entirely 
depends on improvement in the determination of these hadronic factors. 

The quantities of our interest are T^^^ (0) and fs^ \J Be, .The reader is 
referred to the papers [3S] - [iO] and the papers gTJ for lattice and sum 
rule determination of T-^^^ (0), respectively. The corresponding values are in 

the range 0.2 — 0.4. The relevant references for Je^ \J Bb^ are given by papers 
[43] - [46|. Looking at the data one can see that there is no clear agreement 



^"Notice that we have included the _B-meson mass ms to both factors in l|20| l to provide 
normalization. Indeed, it can be understood as being taken from real experiment, on the other 
hand since any reliable lattice simulation must be correctly normalized to this experimental 
value of the mass, the treatment of rriB as a lattice output should make in fact no difference. 
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between, e.g. lattice computations and light-cone sum rules results. Moreover, 
the errors given by the authors of the cited lattice papers are mostly statistical 
ones. The procedure of correct treatment of systematic errors in this case is not 
yet known. In fact, the same is true for the sum rule calculations. In general, 
precise determination of T^^^ (0) on the lattice is very difScult and reliability 
of the calculations done so far is debatable (see [37j and recent discussion in 
[48]). However the utmost importance of this measurement, which hopefully 
will be done in the nearest future on new improved lattices in unquenched case 
cannot be overestimated. 

Thus, having no better strategy at the moment, we will be conservative in 

our error treatment and take for the input value of Jb, \J Bb^ 

Ib, \J~^, = (280 ± 40) MeV 

while we also consider three sets of possible values for Ti~*^ (0), where the 
errors correspond to those reported in the cited papers: 

Set A : T^r^' (0) ^ (0.25 ± 0.05) 
Set B : T^b'^' (0) = (0.30 ± 0.05) 
Set C : T^cT^' (0) = (0.35 ± 0.05) 

The ultimate goal should be to reach the accuracy of the lattice computations 
comparable to the accuracy of the r.h.s. of (|25l) . 

Finally, let us recall the experimental data for Bg meson levels splitting AA/^ 
[26] and branching ratios for the decay B K*j. They are given by 

AM, = [17.3Sto fi{stat) ± 0.07{sys)]ps'^ (22) 

and 

Sr(S" ^ if*"7) ^ (4.25 ±0.31 ±0.24) • 10"^ [49] (23) 

Sr(S- ^ if *-7) = (3.87 ±0.28 ±0.26) • 10"^ [50] (24) 

For the life time of meson we use the value t — (1.652±0.014) ps [51], while 
the masses (in MeV) are given by [5T] 

ms = (5279.0 ±0.5) ; ms^ = (5367.5 ± 1.8) ; rriif. = (891.66 ± 0.26) 

Other short-distance inputs are collected in the Table 1. 

We have all input data now to estimate the ratios Cs,K' ■ According to the 
three choices of numerical value for the form-factor T-^^^ (0) we get three sets 
of Ci The results are presented in Table 2. For graphical presentation one 

can use planar ternary coordinates where the constraint X]i=i I'^sCl*]^* = is 
satisfied automatically. Each solution is represented by a single point on this 
plane with the distance from the i-th axis to the point given by \ogC[\*- It is 
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taken positive for two axes forming an angle the point belongs to and negative 
for the remaining distant axis. With this rule each point on the plane satisfies 
the constraint (fTB|) . Some sample result for the case [44]-B is shown on Figure 
4. Notice, that the bars correspond to la deviation in Ci^lf * > '^ot in log Ci*i<-» • 
The fact that they cross the corresponding axes means less than la deviation of 
the actual result from the SM prediction. The origin of this ternary coordinate 
system corresponds to logCs'^. ^ ^ *' which is the SM solution. 

The main qualitative conclusion is perhaps not surprising: with the rea- 
sonable choice of parameters we observe no evidence for NP within error bars. 
There are two optimistic remarks however. The first is that our errors are very 
conservative and significant reduction of at least some of them is foreseen in 
the nearest future. Secondly, the errors in the Table 2 are not independent. 
There are two sorts of correlations. The first is the uninteresting "kinematical" 
one, following from the constraint (jl6p . The second pattern corresponds to the 

error correlation for lattice simulations of T^^^ (0) and fs, \J Bb^ ■ So far 
these two inputs have been measured independently, by different lattice groups 
and within different procedures. Correspondingly, the errors shown in the Ta- 
ble 2 are also treated as independent. On the other hand, it is reasonable to 
expect an error reduction for the simultaneous calculation of T^^^ (0) and 

Ibs \J ^^^d we call the attention to importance of such simulation, using the 
same framework (lattice action, chiral extrapolation procedure etc) and uniform 
error treatment. It is reasonable to expect that this would result in a better 
accuracy, first of all for the quantity Cl^l-* ■ Speaking differently, if one assumes 
no NP (i.e. Ci — 1) ^o get 

m,\Tf^pp)\ ^ j^^^^ ^ ^,^^ji/2 (25) 

where the uncertainty in the numerical factor 778 is of order 5% and is mostly 
theoretical This SM prediction demonstrates the level of precision the lattice 
computations must reach in order to make reliable conclusion about NP based 
on the lattice results. We consider the check of ([25| on the lattice as a task of 
primary importance. 

3.2 Analysis of the angle a 

The angle a can be extracted from the two-body decay modes of B into light 
hadrons tt, p and K (see recent review [S5]). From the theoretical point of view 
the best channel seems at present to be i? ^ pp [53l [54] . The most promising 
channel for a at LHCb however is B ^ pvr — > tttttt [55l [56] . The basic idea of the 
analysis [57] is to study the interference of the tree amplitude proportional to 
the weak phase factor e'^ from V*ijVud and the penguin amplitude proportional 

^^The uncertainty in experimental value of AMs is small. 
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to factor e~'^^ from V^lVtd- Writing down also the amplitudes for CP-conjugated 
modes and imposing isospin relations, one can fit four amplitudes, four strong 
phases and one weak phase from 11 observables (see details in [Snj). The ex- 
pected uncertainty in a of LHCb is about 10° in one year of running [52]. It 
can be mentioned that the recent result presented by BaBar collaboration [55] 
for a from B —> pir channel is 

a =(114 ±39)° (26) 

while the data uncertainty for B ^ pp mode is ±13° [49 . The above analysis 
assumes no electroweak penguin contributions. According to the estimates j54| . 
Saswp — —1.5°. The isospin breaking effects controlled by parameter (to^ — 
mu)/^QCD are expected to be of the same order of magnitude. 

For a defined as an argument of the amplitude ratio one gets (see details in, 

e.g. m) 



2aeff = arg 



,o,, AiB^f) 
A{B ^ /) 



arg 



-17 _ ie+iSa 
gi7 j.QiO—iSa 



(27) 



where 9i2 is the Bq — Bq mixing angle, 9 is strong penguin phase, r is an 
absolute value of penguin-to-tree ratio and Sa is possible weak NP penguin 
phase In the absence of penguins, i.e. if r = and if 6*12 = 2(3 (as in the 
SM), one gets cte// = atree with atree defined by (|14p . It is worth stressing 
(see early discussion of related issue in [12]) that for the discussed scenario the 
corresponding NP phase shift 6a is to coincide up to a sign with that to the 
angles (3 and x- 

SPnP = P{B^4>Ks) - P{B^J/4,Ks) = = 6x (28) 

due to the assumed J7-spin invarianceF^ Also it has to be noticed that the box 
diagram corresponding to the — B'^ mixing contributes identically to the 
discussed decay modes and its contribution to the phase (with a possible NP 
part) is canceled in (|28|) . Certainly beyond the SM one could have 612 ^ 2/3, 
but this phase shift may have no direct relation to the discussed shift Sa = 
—(5/3, resulted from the penguin process. Thus we are left with the only NP 
contribution from the penguin-mediated decay (with respect to the tree level 
one). We see that the ability to extract 6a from experiment (i.e. from aeff) 
crucially depends on the value of r, since given experimental uncertainty in aeff 
corresponds to larger uncertainty in da smaller the ratio r is. The combined 
fit of the data for B ^ pir and other modes (notably Bq K^po) taking into 
account nonzero penguin NP phase 5a is being performed and will be reported 
elsewhere. Here we would like to notice that the experimental accuracy of 6f3 
is currently limited by the statistics of i? ^ 't'Ks decay and recent update for 
sin 2/3 from penguin decay modes as given by [23; is sin 2/3|peng = 0.58^p pg±0.13 



^^The above assignment is self-consistent for r <J 1, where terms in r which are nonUnear 
in penguin amplitudes can be neglected. 

^''See discussion of the related issue in supersymmetric context in I60| . 
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which correspond to about 15° uncertainty window in the angle (3. It is worth 
mentioning that this penguin-dominated mode would not allow to get sin 2/3 
(and hence (5/3) with competitive precision at LHCb, since the latter is expected 
to be about 0.2 in 2fb"^ of running [22 . Higher accuracy should be possible for 
Super-B factories. 

4 Conclusion 

The standard approach to study CKM matrix is to overconstrain the UT using 
all available experimental information. However not all constrains on the (p, 77) 
plane are sensitive to NP, at least if the latter is taken in the form of next-to- 
minimal flavor violation. Some (such as A/W^/A/Ws) do not distinguish the SM 
from many NP scenarios just by construction, while others (such as relation 
([5])) are insensitive to NP because of the specific profile of the UT (a close to 
90°). In this sense there are two possible points of view regarding the fact that 
up to now all constraints on (p, rf) plane agree with each other. The first one 
is that there are no sizeable NP effects seen in flavor physics. The second one 
is that UT is simply not suitable for the purpose (since NP is not present in 
the angles determined from the tree processes and could also cancel from the 
sides) and the room for manifestations of NP in b-physics observables is in fact 
not so small (because the uncertainties are still rather large). Following the 
latter attitude, we have discussed in this paper a complementary analysis of the 
data on the CKM matrix elements. Its key feature is the use of CKM matrix 
elements ratios which are sensitive to NP provided it contributes differently 
to the processes of different topology. In this sense the quantities Q are 
different from the ratios like AM^/ AMg since the short-distance part is kept 
in the former. Moreover, since we have more than one choice for observables 
a given CKM matrix element is extracted from, we could have relations of 
the form (fTC)) . leaving unconstrained more than one degree of freedom. Thus 
the lattice simulations must match several hadronic inputs simultaneously (and 
not just one). This, we believe, will allow to reduce the corresponding errors 
and consequently to make the proposed probes more sensitive to the NP0 
Speaking differently, one of our main messages to the lattice community is that 
the importance of further reducing uncertainties in the ratio f is limited with 
respect to the calculation of hadronic inputs entering the definition of C's since 
the latter are more sensitive to NP than the former. 

Concerning the determination of UT angles which are free from lattice uncer- 
tainties, we advocate the importance of estimates of the angle Sa corresponding 
to the penguin amplitude extracted from B pir and other modes (and hence 
subject of possible NP shifts). The accuracy of such a comparison can be com- 
parable or better at LHCb than for sin 2/3 extracted from B J/tpKs and 
B —> 4>Ks modes, while the physical meaning is the same; any discrepancy 
between these values would undoubtedly indicate NP. 

^*The importance of correlating of AF = 1 and A_F = 2 processes is stressed in another 
respect in I10| . 
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In principle, nothing prevents one to include the discussed quantities Q^y 
and 5(a,fj,x) into the global fit of the CKM matrix. It is clear that one gets 
essentially no new information in this way, since we deal with the same experi- 
mental observables the standard fitting procedure does. We feel, however, that 
careful analysis of the proposed observables provides an alternative and trans- 
parent way of looking at NP effects. This strategy can become useful in the 
nearest future when LHC data will improve the accuracy of our knowledge of 
the CKM matrix elements dramatically. 
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Table 1: Short-distance quantities from the definition of k. 



Quantity 


Mean 


Error 


Reference 


VB 


0.551 


0.008 


33J 


Ftt 


2.35 


0.06 




rhtimt), GeV 


164.7 


2.8 


EH 


fhb{pole), GeV 


4.65 


0.10 


28J 




0.16 


0.01 


[28] 


mw, GeV 


80.40 


0.03 


f5T] 



Table 2: Numerical results for C^*]^,/AC^*]f* • The abbreviation [43]-A corre- 
sponds to the branching ratio for B — > K*'-f from [43] and the Set A choice for 
T-P^^ (0) = 0.25 ± 0.05, and analogously for other columns. 





[43]-A 


[43]-B 


[43]-C 


[44]-A 


[44]-B 


[44]-C 




0.81/0.21 


0.97/0.23 


1.13/0.25 


0.85/0.22 


1.02/0.24 


1.18/0.26 


A'2) 
^s.K* 


1.12/0.24 


0.94/0.18 


0.80/0.13 


1.07/0.23 


0.89/0.17 


0.77/0.12 


A-i) 


1.10/0.17 


1.10/0.17 


1.10/0.17 


1.10/0.17 


1.10/0.17 


1.10/0.17 
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Figure 1: Error propagation corresponding to the second expression from 
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Figure 2: The same as Fig.l for different values of the angle a. 
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cr(r), ° o-(R,) = 0.02 -0.04 

R,=0.41 




-0.2 -0.1 0.1 0.2 



Figure 3: Error propagation for Rt from ([8]). The curves correspond to a 
0.02 and 0.04. The choice for other parameters is Rb = 0.41, 7=1 rad. 
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Figure 4: The results for logQ j^,, the case [44]-B plotted as a point in ternary 
coordinates. The SM solution is the point at the origin. The algebraic distance 
from the i-th axis is given by logC^']^», positive for two axes forming an angle 
a point belongs to and negative for the remaining distant axis. With this rule 
each point on the plain satisfies the constraint (|16p . The bars correspond to la 
deviation in Ci*l-« ? iiot in log Ci']<-. • The marks on axes set the scale and serve 
mainly for guiding eyes. 
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